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q ' ABSTRACT 

Aims. We re-examine the scenario of X-ray emission from magnetically confined/channelled wind shocks (MCWS) for Ap/Bp stars, 
^"v a model originally developed to explain the ROSAT detection of the AOp star IQ Aur. 

Methods. We present new X-ray observations of the AOp stars a 2 CVn (Chandra) and IQ Aur (XMM - Newton) and discuss our 
findings in the context of X-ray generating mechanisms of magnetic, chemically peculiar, intermediate mass stars. 
Results. The X-ray luminosities of IQ Aur with log L x = 29.6 erg s _I and a 1 CVn with log L x ;£ 26.0 erg s _1 differ by at least three 
orders of magnitude, despite both being AOp stars. By studying a sample of comparison stars, we find that X-ray emission is preferably 
generated by more massive objects such as IQ Aur. Besides a strong, cool plasma component, significant amounts of hot (> 10 MK) 
plasma are present during the quasi-quiescent phase of IQ Aur; further, diagnostics of the UV sensitive f/i line ratio in He-like O vn 
triplet point to X-ray emitting regions well above the stellar surface of IQ Aur. In addition we detect a large flare from IQ Aur with 
temperatures up to ~ 100 MK and a peak X-ray luminosity of log Lx ~ 3 1 .5 erg s . The flare, showing a fast rise and e-folding decay 
time of less than half an hour, originates from a rather compact structure and is accompanied by a significant metallicity increase. The 
X-ray properties of IQ Aur cannot be described by wind shocks only and require the presence of magnetic reconnection. This is most 
evident in the, to our knowledge, first X-ray flare reported from an AOp star. 
*-^Hi Conclusions. Our study indicates that the occurrence the of X-ray emission in AOp stars generated by magnetically channelled wind 

shocks depends on stellar properties such as luminosity which promote a high mass loss rate, whereas magnetic field configuration 
j.^ ' and transient phenomena refine their appearance. While we cannot rule out unknown close companions, the X-ray emission from 

IQ Aur can be described consistently in the MCWS scenario, in which the very strong magnetic confinement of the stellar wind has 
led to the build up of a rigidly rotating disk around the star, where magnetic reconnection and centrifugal breakout events occur. 
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OO ' 1. Introduction hidden in the vicinity of their optically bright primaries have 
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been often proposed as the true origin of the X-ray emission. 



, X-ray emission is found virtually everywhere among the popu- Unresolved companions are a likely explanation for the X-ray 

■ lation of the main-sequence stars. The X-ray generating mech- emission from 'normal' main-sequence stars at these spectral 
; anisms differ however for 'hot' O and early-B stars, and types a nd in some cases they even can be shown to be responsi- 

O ■ 'cool', late-A to M-type stars. In hot stars the correlation be- ble for me X-ray emission, e.g. in th e eclipsing binary systems 

■ tween X-ray and bolom etric luminosity at logL x /L bo i ~ -7 like g OB (ISchmitt & Kursterfl993l) . However, especially pecu- 
t^H , (|Berghoeferet_alJU997|) is explained by X-ray emission aris- liar Ap/Bp stars and young He rbig AeBe s tars are prime candi- 

^ ■ ing from shocks in instabilities in their radiatively driven winds, dates for being intrinsic X -ray emitters (e.g. lBabel & Montmerld 

k> ' In 0001 StalS ' X " lay emission ori g inates in magnetic structures 19971; IStelzer etal]|2009h . In the here discussed Ap/Bp stars an 

^ , generated by dynamo processes, powered by an interplay of interplay between stellar winds and large scale magnetic fields 

H ■ convective motions in the outer layers and stellar rotation. The is thought to play the major role in the X-ray generation. 
- . .' observed activ ity levels range from logLx/^boi ~ -7 ... - 3 
(ISchmittlll997h . whereas an activity-rotation relation is present 

unt il the dynam o saturates (for an overview on stellar coronae Ap/Bp stars are magnetic, chemically peculiar (CP) stars that 
see[Gude] 2004). In intermediate mass stars, around spectral type have completed an often significant fraction of their main se- 
late-A, the outer convection zone becomes increasingly thinner quence lifetime, i.e. they are not especially young. Only a small 
and the dynamo efficiency declines strongly. A prominent ex- fraction of about 5 % of the mid-B to mid-A type stars are mag- 
ample is the A7 star Altair, o ne of the hottest magnetica lly ac- netic; these stars typically exhibit various chemical peculiarities 
tive stars detected in X-rays (Robrade & Schmitt 2009). Thus and a rather slow rotation compared to non-magnetic stars of 
one expects a virtually X-ray dark population at spectral types similar spectral type. The origin of their magnetic field is de- 
mid-B to mid-A, where stars neither possess an outer convec- bated, but often thought to be of fossil origin, consistent with 
tion zone nor drive strong stellar winds. The overall X-ray detec- the finding that only a small fraction is magnetic and that field 
tion rate among those stars is with ~ 10 % indeed quite low, but strength and rotation are independent. The magnetic field geom- 
surpri singly X-ray emission was detected from several of these etry in Ap/Bp stars is dominated by rather simple, large scale 
stars (Schroder & Schmitt 2007). Consequently, low-mass stars structures like a dipole and hence it is fundamentally different 
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from the complex field geomet ry of magneticall y active late-type 
stars (see e.g. the overview bv lLandstreetlll992l) . 

The prototypical case for a presumably single, magnetic 
Ap star with detected X-ray emission is IQ Aur. In a 10 ks 
pointed observation with ROSAT PSPC, IQ Aur was shown to 
be an X-ray source with a high X-ray luminosity of Lx = 
4 x 10 29 ergs -1 and a quit e soft spectrum with Tx = 0.3 keV 
dBabel & Montmerld [l997l BM97). Thus a hypothetical late- 
type companion was a rather unlikely explanation. In order to 
explain the X-ray emission from IQ Aur and, more generally of 
Ap/Bp stars, BM97 introduced the 'magnetically confined wind- 
shock model' (MCWS). In the MCWS model the radiatively 
driven wind components from both hemispheres of the star are 
magnetically confined and forced to collide in the vicinity of 
its equatorial plane. As a consequence, strong shocks form and 
give rise to very efficient X-ray production from kinetic wind 
energy. Although some fine-tuning is required to obtain a rather 
high mass loss rate, the X-ray emission of IQ Aur can actually 
be accounted for with reasonable assumptions on the wind pa- 
rameters; specifically BM97 assumed a stellar wind with V M 
800 km s and a mass loss rate of M Q = 10~ 10 - 1O~ U M yr _1 . 

Notable advancements of the static BM97 model were 
achieved by magnetohydrodynamic simul ations of magneti- 
cally channelled line-driven stellar winds (lud-Doula & OwockH 
2002). In a further study, these MHD simulations investigate 
an extended p arameter space and include stellar rotation in 
the modeling dud -Doula et alj 120081) . Dynamic variants of the 
MCWS model were successfully applied to Herbig AeBe and O- 
type stars and may actually extend to a whole class of magnetic 
early-type stars. In this respect it can be considered as a 'stan- 
dard' model that is invoked to explain a variety of phenomena 
like X-ray overluminosity, hard spectral components, flares or 
rotational modulation. However, re cent X-ray st u dies o f Ap/Bp 
stars often remained inconclusive. IStelzer et al.l (120061) studied 
X-ray detected intermediate mass main-sequence stars (B5-A0) 
with Chandra including three CP stars, but they do not find 
strong support for the MCWS model, especia lly due to the non- 
detect ion of the magnetic Ap star HD 13388Q. ICzesla & Schmittl 
J2007h studied mid-B to mid- A magnetic stars; they detected 
X-ray emission from the AOp star GL Lac ( Babcock's star), 
though with a much higher X-ray temperature and luminosity 
than IQ Aur, but also report the non-detection of the AOp star 
HD 184905. They find that the presence of a magnetic field is 
not a sufficient criterion for X-ray emission; further X-ray lumi- 
nosity and field strength are not correlated. 

These findings make it desirable to investigate the MCWS 
scenario in its 'original' environment and we present contempo- 
rary X-ray observations of two prototypical AOp stars, a 2 CVn 
and IQ Aur. Our paper is structured as follows: we summarize 
the targets properties in Sect.|2] in Sect. [3] we describe observa- 
tions and data analysis, present our findings in Sect. |4] and dis- 
cuss them in a broader contexts in Sectj5] 



Table 1. Basic properties of our targets. 



Target IQ Aur 

X-ray obs. XMM 

Dist. (pc) 137 

Sp.(chem) type AOp - CP2, Si 

log L(L ), log T(K) 2.5 ±0.1, 4.16: 

M (M Q ), R (R e ) 4.0 ±0.1, 2.9 ±0.2 

log t (yr), t (frac.) 7.8 (7.5-7.9), 0.4 

v sini (km/s), P m (d) 56, 2.47 

S cft (kG) 0.74 ± 0.4 



a 2 CVn 

Chandra-ACIS-S 
34 

A0p-CP2, SiEuHgCr 
0.01 2.0 ±0.1, 4.06 ±0.02 
3.0 ±0.1, 2.5 ±0.3 
8.3 (8.1-8.4), 0.52 
18, 5.47 
1.35 ± 0.4 



find that the overall magnetic structure of a 2 CVn can be mod- 
elled with an 4.6 kG dipolar field, an inclination of i « 120° 
and an obliquity, i.e. the angle between the rotational and the 
magnetic axes, of ft ~ 80°. Nevertheless, also higher complexity 
in the magnetic topology is present. X-ray emission from this 

-_±.0 keV) at a level of 

28.8 ergs" 1 



( Schmitt et al.l 1985 


) and ROSAT (0.1- 


gL x =28.6 ergs" 1 ( 


Schroder & Schmitt 



logix 

2.4 ke V) at a level of log L x 
2007), but it was not resolved in both observations and the con- 
tribution from the individual components remained unknown. 

The field AOp star IQ Aur (HD 34452, HR 1732) belongs to 
the class of a 2 CVn spectrum variable stars and is thought to be 
a single star with V=5.4 and small £b-v — 0.01. BM97, who re- 
port the ROSATX-my detection at log Lx = 29.6 erg s _1 , assume 
a 4 kG dipolar field for IQ Aur and determine a n inclination 
of i » 32° and an obliquity of /3 « 50 °, whereas Hubrig et al.l 
d2007b find i = 87° and of B = 7°. iBohlender et all i 19931) 
state: 'a large range of values for i is possible: if i is large ft 
is fairly small and if i is as small as allowed by our assump- 
tions, (3 may be more than 80°. Also stellar parameter somewhat 
differ, e.g. BM97 assume log L/L = 2.9 ± 0.2, M _ 4 8 M» 
R = 5.1 + 1.1 and V rot = 105 kms 1 whereas lHubrig et al.l(l2007l) 
give log L/L Q =2.3 ±0.1, M = 3.2 ± 0.2 M , R = 2.2 + 0.4 R Q 
and V mt = 46 kms~'. Since the adopted distance is the same, the 
stellar data and viewing geometry of IQ Aur appear to be only 
moderately constrained and ti me- variability might p l ay an addi- 
tional role. The values from lKochukhov & Bagnula (12006). that 
are commonly used throughout the paper, are intermediate ones. 

B oth AOp stars are CP2 stars in the classification scheme 
from iPrestonl (1 19741) and might be expected to be overall simi- 
lar given their identical spectral type. However, IQ Aur is more 
luminous, more massive and younger than a 2 CVn. Actually, 
IQ Aur is extreme for an AOp star; as a main-sequence star 
it would be with its B — V = -0.16 of spectral type mid- 
B. The adopte d basic stellar data for I Q Aur and a 2 CVn a s 
compiled from | Bohlender et al.1 (1 1 993b. iBvchkov et al.l (120031) . 
IKochukhov & Bagnuld d2006l) and lKochukhov & Wadd (2010) 
are summarized in TableQ] 



2. The targets a 2 CVn and IQ Aur 

The AOp star a 2 CVn (HD 112413, HR 4915) is the primary 
component (V =2.9) in a wide binary (angular separation 21.3 ") 
with the F0V star HR 4914 as secondary, a 2 CVn is the pro- 
totype of a class of spectrum and brightness variable stars, 
caused by the interplay of stellar rotation and starspots, i.e. 
metal inhomogeneities in the stellar atmosphere. These stars 
have spectral classifications from late-B to mid-A type and 
belong to the class of ch emically peculiar, magnetic stars. 
IKoch ukhov & Wade] d2010h study spectropolarimetric data and 



3. Observations and data analysis 

We observed a 2 CVn in March 2009 with Chandra ACIS- 
S for 15 ks (Obs.ID 9923). The AOp star IQ Aur was ob- 
served in February 2010 with XMM-Newton for 1 15 ks (Obs.ID 
0600320 1 1 ) in the full frame mode with the thick filter inserted. 
Due to the optical brightness of the target the OM (Optical 
Monitor) had to be closed. 

We utilize the software packages CIA04.2 and SAS 10 to 
process and analyze the data, apply standard selection criteria 
and use the rather well calibrated energy ranges 0.2-8.0 keV 
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Fig. 1. Left panel: ACIS-S image of the a 2 CVn, the AOp star 
HR 4915 is X-ray dark. Right panel: MOS1 image of the AOp 
star IQ Aur. The circles denote the respective optical positions. 

(XMM-Newton) and 0.3-8.0 keV (Chandra). Concerning the 
XMM-Newton observation of IQ Aur, For the spectral analysis 
of the XMM-Newton EPIC (European Photon Imaging Camera) 
data of IQ Aur source photons are extracted from a 30" circular 
region around the source, the background is taken from nearby 
source-free regions on the MOS and the PN detectors respec- 
tively. Of the EPIC detectors the PN provides a higher sensi- 
tivity, whereas the MOS provides a slightly better spatial and 
spectral resolution. The signal in the RGS (Reflection Grating 
Spectrometer) detectors is rather weak, we therefore extracted 
high resolution RGS spectra from a 90 % PSF source region and 
performed a cross-check with a narrower 70 % PSF extraction 
area, in which the background is further suppressed. An addi- 
tional problem are effective area drops in the RGS that affect e.g. 
the region around the O vn intercombination line where several 
so-called 'bad' channels lead to a reduced number of detected 
counts. To measure the line flux es we use the fitting program 
CORA dNess & Wichmannl2b 02). that accounts for these detec- 
tor inhomogeneities. 

We perform spectral anal ysis of the EPIC/ ACIS-S data with 
XSPEC V12.3 (lArnaudl[T996l) . To model the spectra we use 
photo-absorbed, multi-temperature APEC models with abun- 
dances relative to solar values as given by iGrevesse & Sauvall 
(1998), errors from spectral modeling denote 90 % confidence 
level. We note that the interdependence of absorption, metallicity 
and emission measure (EM = J n e nj]dV) may result in different 
models of similar fit quality with the consequence that absolute 
abundances are more poorly constrained than abundance ratios 
or relative changes. Also the incomplete knowledge of atomic 
data used to calculate the spectral models adds to these uncer- 
tainties, however they do not alter the main conclusions of this 
work. 

4. Results 

Here we report on the results relevant to the origin of X-ray emis- 
sion in Ap/Bp stars and present an analysis of the IQ Aur data, 
model its spectra, take a detailed look at the flare and investigate 
X-ray lines from the high-resolution spectrum. 

4.1. Source detections, light curves and a flare 

The ACIS-S image of the a 2 CVn field, shown in Fig.Q] clearly 
reveals that the AOp star a 2 CVn (HR 4915) is X-ray dark. In 
the 0.3-3.0 keV band we detect no photon in a 2 "source region 
(95 % PSF) centered on the expected optical position at an esti- 
mated background of 0.2 photons. We detect the F0 star HR 49 14 
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Fig. 2. X-ray light curve of IQ Aur from the summed EPIC data 
in the 0.2 - 8.0 keV band. Inset: The flare on a logarithmic scale 
in three energy bands as seen by the PN detector. 
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Fig. 3. Quasi-quiescent part of the X-ray light curve, phased with 
the rotation period of IQ Aur and binned to 1 ks (asterisks) and 
3 ks (diamonds). The average rate before and after the flare is 
indicated by a red dashed line respectively. 



with an X-ray luminosity of log Lx = 28.6 erg s _1 , a rather typi- 
cal value. For its position we find a small offset of about 0.5 " be- 
tween the X-ray and expected optical position; nevertheless, ap- 
plying a spatial correction to HR 4915 does not change the out- 
come. Adopting Poissonian statistics we derive an 95 % conf. 
upper limit of three counts for the source region. Assuming a 
plasma temperature of 0.3 (0.8) keV this corresponds to an up- 
per limit on the X-ray luminosity of Lx < 0.6 (1 .0) x 10 26 erg s" 1 , 
three to four orders of magnitude below the X-ray luminosity of 
IQ Aur (Lx — 4 x 10 29 ergs -1 ), despite both being apparently 
similar stars. The flux limit for a 2 CVn converts into the by far 
the lowest upper limit obtained on the X-ray brightness for an 
AOp star. 

In contrast, IQ Aur is detected in all XMM-Newton detec- 
tors; see the right panel of Fig.Q] where we show the MOS1 
image. The derived X-ray position is about < 1 " off from the 
expected optical position, a value within the positional uncer- 
tainty for XMM-Newton. To investigate the X-ray variability of 
IQ Aur, we extracted photons from a 35 "region (~ 90% en- 
circled energy) around the X-ray source position on the EPIC 
detectors. 

In Fig. [2] we show the background subtracted light curve 
with 100 s time binning. Beside minor variability that is present 
throughout the observation, the X-ray count rate increases after 
49 ks by a factor of roughly 40 within less than half an hour, 
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the total flare event has a duration of about two hours. In the in- 
set of Fig. |2] we show the light curve of the flare part as seen by 
the PN detector, also separately in a soft and hard energy band. 
The shape of the light curve with its fast rise time of 1000 s 
and e-folding decay time of ~ 1200- 1400 s points to a single, 
rather compact flaring structure. Since for X-ray emitting struc- 
tures/loops the radiative energy loss scales with the square of the 
particle density, compact, i.e. small and dense structures, radi- 
ate their energy much faster away than larger tenuous structures 
of comparable X-ray luminosity. Further, the harder X-ray emis- 
sion increases by a larger factor and its peak precedes the one 
of the softer emission, a behavior that is reminiscent of coro- 
nal flares. We do not detect any spatial shift between the pho- 
tons detected during the quasi-quiescent and the flaring phase 
at an accuracy of < 0.5". In Fig.[3]we show the quasi-quiescent 
part with two different time binnings. Notably, the average quasi- 
quiescent X-ray count rate has increased by roughly 20 % higher 
after the flare and the light curve shows slightly stronger vari- 
ability, although the spectral hardness given by the ratio of the 
count rates in the 0.2-0.8 keV vs. 0.8-2.0 keV energy band is 
virtually identical compared to the pre-flare phase. 

We do not see a strong modulation of the X-ray light curve 
that might be expected in an oblique rotator model (P rot = 
213 ks). Significant, roughly sinusoidal shaped rotational modu- 
lation would be present, if the star partially occults X-ray emit- 
ting material depending on rotational phase or when a disk 
around the star is viewed under different angles. An example 
are the X-ra y brightness variat ions of the massive magnetic 
star l OriC (Gag ni"e"t al.ll2.Q05h . This model assumes cylindri- 
cal symmetry for the system as is likely appropriate for a star 
with a disk and what is commonly observed when e.g. mag- 
netic field variations are monitored, but clearly we cannot com- 
pletely exclude unexpected phenomena during the not covered 
phase. Brightness variations on timescales of 20 ks (~ 0. 1 phase) 
between consecutive bins during the pre-flare or the post-flare 
phase have an amplitude of < 10 % and appear rather irregular. 
However, the strength of the modulation depends on the plasma 
location and the geometry of the system and for a small obliquity 
the effects are likely minor in any case. 

4.2. Spectral properties in quasi-quiescence 

The global spectral properties of IQ Aur are determined from the 
EPIC spectra by using photo-absorbed multi-temperature mod- 
els with free metallicity. Given the large variations in count rate 
over the observation, we separate the data into a flaring phase 
(PN: 49.1-58.1 ks) and a quasi-quiescent (QQ) phase for spec- 
tra analysis. As an example we show the respective PN spectra 
with the applied spectral models in Fig.|4j the derived spectral 
properties from our modeling are summarized in Table |2] 

The description of the quasi-quiescent spectrum clearly re- 
quires, beside the 0.3 keV component seen by ROSAT, the pres- 
ence of additional hotter plasma components. Applying a one 
temperature model for example to the MOS spectra results in an 
average temperature of kT = 0.75 keV, albeit with > 2 
the fit is quite poor. Models with three temperature compo- 
nents describe the individual data sets rather well. The absorp- 
tion component required to model the X-ray spectra of IQ Aur 
is rather weak, compatible with the findings from the ROSAT 
data. We again find strong emission from rather cool plasma 
at ~ 3 MK, but overall the plasma temperatures derived from 
the XMM-Newton spectra are significantly higher than those de- 
termined with ROSAT, where only weak indications for hot- 
ter plasma were found due to the limited sensitivity of the in- 
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Fig. 4. XMM/PN spectra of the quasi-quiescent (black) and the 
flaring (red) phase with applied spectral models. 



strument. We detect a significant contribution from > 10 MK 
plasma even in the quasi-quiescent phase of IQ Aur, here it ac- 
counts for about 30 % of the total emission measure. We find 
that a subsolar metallicity best describes the data. Using so- 
lar abundances reduces the emission measure, but also worsens 
the model quality (x 2 red = 1-24), albeit the temperature struc- 
ture remains nearly identical. Allowing individual abundances 
to vary independently does not lead to robust results, therefore 
a global metallicity is used in all models. The best fit model 
corresponds to a source flux of 1.7 x 10~ 13 ergcm^s -1 and 
an X-ray luminosity of Lx - 4.2 x 10 29 ergs -1 in the 0.2- 
2.0 keV band; for the ROSAT 0.1-2.4 keV band we derive a 
roughly 20 % higher flux. This corresponds to an activity level 
of logLx/^boi ~ -6.4, assuming that the X-rays are emitted by 
IQ Aur. Considering the ROSAT detection at a very similar X-ray 
luminosity of Lx = 4.0 x 10 29 erg s _1 nearly 15 years ago, IQ Am- 
is likely emitting X-rays at a rather constant level for more than 
a decade. 

4.3. Global flare properties 

During the flare phase, IQ Aur has not only significantly bright- 
ened, but as shown in Fig.|H also the associated spectrum is 
much harder. The flare is most prominent at high temperatures, 
but large excess emission is also present at lower temperatures. 
During the peak of the strong flare the emitted X-ray luminosity 
of IQ Aur rises to about L x = 3.2 x 10 31 erg s _1 (0.2- 10.0 keV), 
corresponding to a flux increase of nearly a factor 100 compared 
to the quasi-quiescent state. With an average excess luminosity 
of L x = 6.2 x 10 30 ergs -1 and a duration of 9000 s, the flare 
released about Ex = 5.6 x 10 34 erg only at soft X-ray energies. 
The flare emission is dominated by extremely hot plasma and 
even the average temperature is found to be around 45 MK. In 
the PN spectrum the prominent emission line complex from He- 
like Fexxv is clearly visible at 6.7 keV, consistent with a hot 
plasma component at temperatures of around 60-70 MK as de- 
duced from the spectral model. Weak excess emission is present 
at the low energy tail of this feature, that could be due to fluores- 
cent emission from neutral iron from the Fe Ka line at 6.4 keV. 
However, formally the line is not detected at the 90 % CL and 
its contribution to the X-ray flux is with an equivalent width of 
EW< 100 eV at best marginal. No significant changes in the 
absorption component are detected; in contrast, the metallicity 
of the flare plasma is found to be significantly higher than those 
of the quasi-quiescent phase. While the absolute scale is rather 
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Table 2. Spectral models of IQ Aur (EPIC data). 



Par. 


QQ 


Flare 


Unit 


N a 


1.1+0.9 


1.0+0.9 


10 20 cnT 2 


Abund. 


0.31+0.08 


0.67+0.28 


solar 


T\ 


0.26±0.03 


0.71±0.09 


keV 


EM\ 


1.75+0.71 


4.7±2.0 


10 52 cnr 3 


T 2 


0.64±0.15 


1.59±0.19 


keV 


EMj 


1.58±0.56 


14.4±6.9 


10 52 cnT 3 


T 3 


1.10±0.11 


6.18±1.46 


keV 


EM 3 


1.88±0.66 


22.4±7.1 


10 52 cnr 3 




1.09 (413) 


0.97 (555) 




L x (em. 0.2- 10.0 keV) 


4.5 x 10 29 


6.7 x 10 30 


erg s _1 



poorly constrained, we find a distinct metallicity trend, most no- 
ticeable when further separating the flare into individual time 
intervals as discussed in Sect l4.4l 

4.4. The flare under scrutiny 

To study the evolution of the flare and its underlying structure in 
greater detail, we perform a time resolved analysis for each of 
the six time segments as shown in the upper panel of Fig. [5] To 
derive the properties of the flare plasma (1-6), we added plasma 
components to the quasi-quiescent (QQ) model derived for the 
PN data. Here we use models with one or two temperature com- 
ponents and with or without metallicity as a free parameter to 
describe the flare plasma. The models lead to overall consistent 
results, see for example the time-evolution of the plasma proper- 
ties presented in the middle panel of Fig.[5]for models with solar 
metallicity. The temperature given for the 2-T model is the emis- 
sion measure weighted temperature. The flare reaches around 
maximum an X-ray brightness of about logLx = 31.5 ergs -1 
and an average plasma temperature of around 80 MK, shortly 
after the impulsive heating during the rise phase (bin 1), even 
roughly 100 MK plasma is present. The maximum temperature 
leads the maximum emission measure and thus in general the 
flare shows the temporal evolution of a typical coronal flare, i.e. 
a clock-wise turn in the temperature vs. density plane (see e.g. 
Idudell 12004 . Here the square-root of the emission measure is 
used as a density proxy, thus a fixed spatial extend of the X- 
ray emission region is assumed. While the initial decay is ap- 
parently steep, but only poorly constrained (bins 2-3), the flare 
evolution trajectory shows some flattening (bins 3-6) of the de- 
cay path. The origin of this behavior is unknown; it could be due 
to minor re-heating events in the same structure or produced by 
additional reconnection events in other magnetic structures, that 
again might be triggered by the large flare. Another possibility 
is a change in size of the emitting region over the duration of the 
event. 

Assuming a single, loop-like emitting structure we can de- 
rive its size by adopting the formalism described in iReale et all 
(Il997l) and p arameters applicable to XMM/EPIC as given in 
Realg (120071) . By using the maximum plasma temperature and 
the flare decay time obtained from the light curve one can de- 
termine the size of the flaring structure, that is treated as an 
instantaneously heated coronal loop. Basically one finds that 
L °c Tic VTmax, L being the loop half-length, r max the max- 
imum plasma temperature and t the decay time from the X-ray 
light curve and an additional scaling factor. This factor depends 
on the amount of additional heating during the decay, thus one 
needs to estimate the importance of sustained heating. This is 
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Fig. 5. Time resolved study of the IQ Aur flare with XMM/PN 
data. Top: time intervals used in the analysis. Middle: tempo- 
ral evolution of the plasma in the temperature/density plane 
(black/red: 1/2-T model). Bottom: metallicity vs. plasma tem- 
perature. 



done by fitting the decay path derived from the one tempera- 
ture model with a straight line to about about one tenth of the 
maximum intensity. This corresponds to bins 2 - 5 in the middle 
panel of Fig. [5] where we derive a slope of ~ 1 for the decay and 
thus moderate re-heating has to be taken into account. Adopting 
the appropriate modeling parameter, i.e. r o b s = 65 MK, an e- 
folding decay time of r = 1200 s (measured to 1/10 of the 
maximum) and a correction factor for a slope of one, we obtain 
L = 1.8 X 10 10 cm (0.1 R») for the loop half-length. 

Since a coronal loop model might not be perfectly appropri- 
ate for IQ Aur and possible heating at later phases of the flare is 
only poorly constrained, we use as a conservative approach our 
most extreme values, i.e. 100 MK plasma temperature and 1400 s 
decay time and ignore ongoing heating. This method results in 
an upper limit for the size and we derive L < 3.5 x 10 10 cm 
(0.2 R*) as maximum extent of the structure. While the flaring 
structure might be much smaller, even the upper limit corre- 
sponds to a moderately sized structure when compared to the 
stellar dimensions of IQ Aur. This finding excludes very large 
magnetic structures, e.g. from the global dipole, as origin of the 
flare. 

The IQ Aur flare is accompanied by a significant change in 
the metallicity of the X-ray emitting plasma as illustrated in the 
lower panel of Fig. [5] where we plot the results obtained from 
a 1-T model with metallicity as free parameter. During the im- 
pulsive phase of the flare as traced by the plasma temperature, 
more metal rich material is heated to X-ray temperatures. The 
freshly heated plasma from the large flare exhibits solar or even 
slightly super-solar metallicity and thus a two to three times 
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Fig. 6. High resolution X-ray spectrum of IQ Aur. Top: Total 
RGS 1 +2 merged and flux-converted. Bottom: O vn triplet with 
model accounting for effective area inhomogeneities. 



higher metallicity than the X-ray emitting plasma observed in 
the quasi-quiescent phase. While the absolute scale is only mod- 
erately constrained as outlined in Sect,[3] the relative changes 
are quite robust and also found in two temperature models or the 
analysis the MOS spectra. In the final stages of the flare evo- 
lution the remaining material has cooled down to temperatures 
below 1 keV, exhibits pre-flare metallicity and becomes rather 
indistinguishable from the quasi-quiescent plasma. Overall, we 
find a clear correlation between plasma temperature and metal- 
licity, indicating a strong modification of the chemical compo- 
sition during the event. Notably, the derived chara cteristics are 
quite typical for a stellar flare (see e.g. Gudel 2003). 

4.5. The high resolution X-ray spectrum 

The high resolution RGS spectrum from IQ Aur has unfortu- 
nately only moderate signal to noise and only rather strong lines 
were clearly detected even in the total spectrum; see Fig.|6] 
where we show a flux-converted spectrum of the merged RGS 
detectors from a 90 % PSF extraction region. The individual line 
fluxes are measured from the total count spectrum and cross- 
checked by using a broad and a narrower 70 % PSF extraction 
region. The lower panel of Fig. [6] shows the O vn triplet with 
applied model, confirming that the forbidden line (22.1 A) is 
significantly stronger than the intercombination line (21.8 A). 
The S/N ratio is quite poor, even in the total spectrum that in- 
cludes the flare. We also investigated the quasi-quiescent and 
flare phase separately, but except stronger emission from 'hotter 
lines' (e.g. Fexvn) during the flare, a quantitative comparison 
suffers from the low S/N. The fluxes of emission lines used in 
the subsequent analysis as determined from the total spectrum 
are given in Table [3] They are corrected for an absorbing col- 
umn of log Nh = 20 cm -2 and we expect a flare contribution of 
roughly 30 % to the given line fluxes. 
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Fig. 7. Show above is the quasi-quiescent O vm/O vn ratio of IQ 
Aur in comparison to accreting/wind driving young stars and 
main sequence stars. 



We utilize the He-like Ovn triplet to determine the dis- 
tance of the X-ray emitting plasma to the surface of IQ Aur. 
Ovn traces the location of plasma at temperatures around 
2 MK and its forbidden (/) to intercombination (i) line ratio 
is sensitive to density and UV radiati on, specifically one find s 
f/j =R /(1+ 4>/<t> c + njn c ), see e.g. lMewe & Schriiverl(ll978l) : 



Pora uet et all (l200ll) . In the low density case (n e < 10 10 cm -3 ) 
additional photons in the /-line are only due to the radiation 
term <p/</> c . The radiation term describes the strength of the UV 
flux at 1630 A and thus depends on the effective temperature 
and distance to the star. Using the fluxes in Table [3] we obtain 
f/i - 3.7 + 3.0, unfortunately the errors are relatively large. The 
ratio is consistent with virtually no radiation field being present 
(Ro = 3.95), but much lower values are formally also allowed. 
Adopting f/i > 1, a photosphere with r e jf = 14450 K and evalu- 
ating the radiation term, one finds that the X-ray emitting region 
is either shielded against the stellar UV-field or the X-rays origi- 
nate predominantly from >1R* above the surface of IQ Aur. As 
an alternative approach, we use the UV flux as measures by the 
IUE satellite and calculate 7^ from the observed spectrum. The 
1630 A UV flux from IQ Aur varies by a factor of nearly three 
over several years; using the lower flux measurements we find 
that the X-ray emitting region has to be located above the star at 
distances of d > 5 R*. This result is independent of our assump- 
tions; if higher densities are present, even larger distances to the 
stellar surface are required. Similarly, applying a correction for 
the weak extinction towards IQ Aur or adopting a higher UV flux 
would again shift the X-ray plasma location to larger distances. 
Therefore we can virtually exclude X-ray emitting structure on 
or close to the star and thus the X-rays are produced at least sev- 
eral stellar radii above the surface of IQ Aur. 

We also investigate the quasi-quiescent and flare phase alone 
and find for both cases a similar large ///-ratio, however with 
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substantial error. Similarly, we study the Ne ix-triplet and find 
an f/i ratio that is slightly higher but again fully consistent with 
the low density/no UV-field limit, however no additional infor- 
mation is gained to O vn. Further, we use the density constraint 
derived from the O vn // /-ratio, as a very rough cross-check on 
the size of the emitting volume. Adopting f/i > 1 implies a den- 
sity of n e < 5 x 10 10 crrT 3 , therefore from this value we can now 
formally derive a minimum size of the emission region. The av- 
erage flare EM is 4 x 10 53 cirT 3 with EM » 0.85 n 2 V and thus 
one obtains for an emitting sphere a radius of r ~ 0.14 R„. We 
again find a rather small structure, in agreement with the estima- 
tions L < 0.2 R„ derived from the above analysis of the flare. 

By using temperature sensitive Ovm/Ovn line ratios, we 
test if additional cool plasma compared to pure coronal sources 
is produced by IQ Aur (see Pig-EJ. The O vn lines form pre- 
dominantly at temperatures around 2 MK, while Ovm traces 
hotter plasma at 3 - 5 MK. Additional O vn emission may be 
produced by accretion or wind shocks, as commonly observed 
in young T Tauri and Herbig AeBe stars, that generate rather 
cool plasma compa red to coronal sources of c omparable X-ray 
brightness, see e.g. iRobrade & Schmittl (120071) . The Ovm/Ovn 
ratio measured for the quasi-quiescent state IQ Aur is slightly 
below the values found for coronal sources, i.e. its X-ray emis- 
sion is only moderately soft for its X-ray luminosity. However, 
it is well above the ratios found for HAeBe stars; young inter- 
mediate mass stars that drive outflows and jets with velocities of 
a few hundred km/s. These stars produce X-ray emission from 
predomin antly cooler plasma, like ly related to shocks in their 
outflows dGunfher & Schmittll2009h . 

5. Discussion 

Our X-ray observations of the AOp stars IQ Aur and a 2 CVn 
present quite a complex X-ray picture. Here we discuss some im- 
plications imposed by our findings in the framework of the mag- 
netically channelled wind shock (MCWS) scenario vs. a possible 
low-mass companion. 

5. 1 . MCWS models for Ap/Bp stars 

To detect MCWS generated X-ray emission from AOp stars, 
some basic conditions have to be fulfilled. The shock-speed has 
to be large enough to heat the plasma to X-ray emitting temper- 
atures, the mass loss rate has be high enough to produce signif- 
icant amounts of hot plasma to be detected with current instru- 
ments and a magnetic field configuration has to be present that 
effectively supports the conversion of kinetic wind energy to X- 
ray emission. 

A key parameter in the dynamic MC WS model is the m ag- 
netic confinement parameter j] t (lud-Doula & Owockill2002l) . It 
basically describes the ratio between the energy density of the 
magnetic field and the kinetic wind energy density and depends 
on stellar properties as given by 77* oc B 2 q R 2 l MVoo', here B eq is 
the equatorial field strength, R t the stellar radius, M the mass 
loss rate and the terminal wind speed. Adopting values for 
IQ Aur, B eq = 2 kG, R, = 2.9 R , M = 10" 11 M Q yr l and 
Voo = 700 km s _1 , gives a confinement parameter in the order of 
77* ~ 10 6 . This value has significant uncertainty, but definitely 
the requirements for strong confinement, i.e. t]* » 1, are ful- 
filled. While at first analytic studies of rigidly rotating magne- 
tospheres indicated that in the case of strong magnetic confine- 
ment a n equatorial disk builds up from corotation radius out- 
wards dTownsend & Owo cki 2005), its formation is also present 



in MHD simulation when including effects of stellar rotation. 
These models predict the formation of a rigidly rotating disk 
by the accumulation of wind material, however episodes of in- 
fall and breakouts limit the build-up of the disk and lead to a 
dynamic but quasi-stationary behavior (ud-Doulaet al . 2008); 
moreover the location of the temporary confined plasma around 
the equatorial plane can be estimated with the formalism pre- 
sented in these works. For moderate rotation and very strong 
confinement the inner edge of disk is roughly given by the 'as- 
sociated Kepler-radius' Rk = W~ 2 ^ with W - V r ot/V i-b and 
V or b = ylGMfRl; using V mt = 50 kms" 1 we obtain R K ~ 5R*. 
The outer edge is determined by the extent of the closed mag- 
netosphere Rc ~ O.IRa, slightly below the Alfven radius given 
by R A ~ y]J A x R t . For IQ Aur we obtain R A w 30 R t \ clearly 
we are in a regime where R A > R^ and thus the conditions re- 
quired by the model for the formation of a rigid disk are fulfilled. 
The disk of IQ Aur is expected to be roughly located between 
~ 5 - 20 R* and breakout events are mainly launched at dis- 
tances of « 20 - 30 R*. We caution that these models assume 
an aligned dipole and were derived for more massive stars, thus 
an extrapolation to our Ap/Bp stars may not be straightforward. 
Nevertheless, the model is generally applicable for magnetically 
channelled line-drive stellar winds and thus appears sufficient for 
the discussion of the overall phenomenology of X-ray emission 
presented in the following. 

5.2. Why is IQ Aur X-ray bright and a 2 CVn X-ray dark? 

The X-ray luminosity produced by MCWS follows the relation 
Lx oc MVoofi*' 4 , whereas the plasma temperature is given by 
T x * 1.15 x 10 5 K(V sh /T00 kms" 1 ) 2 (BM97). Adopting similar 
wind speed and mass loss rate, one expects a 2 CVn to be com- 
parably bright or even moderately X-ray brighter than IQ Aur. 
However, it is at least a factor thousand X-ray fainter than ' ex- 
pected' from the simplified assumption. This implies that either 
the mass loss rates differ by orders of magnitudes or that shock- 
speeds for a 2 CVn are so low, that the plasma does not reach X- 
ray temperatures. Since stellar wind speeds are likely compara- 
ble and strong confinement is achieved for both stars, a very low 
plasma temperature of 7x < 3 x 10 5 K as would be required to 
explain the tight upper limit on Lx for a 2 CVn, is virtually ruled 
out. It also appears unlikely that small differences in the mag- 
netic topology of both stars lead to completely different wind 
channelling and wind shock configurations. A more reasonable 
explanations would be to propose a different mass loss rate and 
that IQ Aur differs by intrinsic properties that promote mass 
loss and thus X-ray generation, for example it is more luminous 
and hotter than a 2 CVn. In this case the X-ray generation would 
depend very sensitively on the underlying stellar parameter(s). 
Alternatively, the MCWS phenomenon could be a transient one 
and currently only IQ Aur is in an active wind phase with a high 
mass loss rate. Unfortunately, mass loss rates of Ap stars usually 
cannot be independently measured, but are determined via mod- 
eling of e.g. their X-ray properties. However, the strong wind 
phase is believed to be connected to abundance anomalies, that 
are also observed in a 2 CVn. It was already noticed by BM97, 
that any abundance anomalies should be removed by the rather 
high mass loss rate, thus requiring that IQ Aur only recently, i.e. 
before a few Myr, entered the active wind phase. Although the 
details of the onset of the strong wind phase are rather unknown, 
they estimate that Ap stars would spends a few percent of their 
lifetime in the strong wind phase, thus explaining that most Ap 
stars are no X-ray sources, at least at the moment. Nevertheless, 
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Table 4. X-ray observations of the comparison Ap/Bp stars. 
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Fig. 9. Correlations for early-Ap/late-Bp stars, Lx vs. stellar 
mass, eff. temperature, eff. magnetic field and age (symbols as 
in Fig.[8}. 



Target 


Sp.type 


Dist. 


Instr./Obs.ID Dur. 


logL x 






(pc) 




(ks) 


(ergs" 1 ) 


HD 12767 


B9.5 Si 


110 


ACIS-S (5391) 


3 


29.4 


HD 27309 


AO SiCr 


96 XMM (0201360201) 


41 


<27.0 


HD 122532 


B9Si 


169 


XMM (0302900101) 


130 


<27.3 


HD 133880 


B9Si 


126 


ACIS-I (2543) 


3 


<28.0 


HD 142301 


B8IIIp 


139 XMM (0142630301) 


20 


28.1 


HD 143473 


B9 Si 


123 


ACIS-S (5385) 


15 


<27.2 


HD 144334 


B8 Si He-w 


149 


ACIS-S (5393) 


3 


30.3 


HD 146001 


B8 He-w 


140 


ACIS-S (5394) 


3 


28.9 


HD 147010 


B9 SiCrFeSr 


143 


ACIS-S (5386) 


15 


<27.3 


HD 184905 


AO SiSrCr 


165 


ACIS-S (5387) 


15 


<27.4 



additional criteria might be necessarily fulfilled to make an AOp 
star a bright X-ray source. 

5.3. X-rays from similar Ap/Bp stars 

Given the very different findings for IQ Aur and a 2 CVn, we 
searched the Chandra and XMM-Newton archives for observa- 
tions of similar Ap/Bp stars, to derive some information on their 
general X-ray properties. As a comparison sample we choose 
late-B to early-A magnetic CP stars with well determined dis- 
tance and stellar parameter, that are either likely single stars or 
spatially resolved in X-rays, although undiscovered binaries in 
the sample are not ruled out. However, their contribution is likely 
minor, since even the overall fraction of close, i.e. spectroscopic, 
binaries in Ap stars of CP2/Si type is w ith ~ 20 % significantl y 
lower than those in non-magnetic stars (Ab t & Snowdenl["l973l) . 
The comparison sample is specified in Tabled some of these ob- 
jects have already b e en presented, albeit with a different focus in 
Istelzer et al.l(l2006l) : ICzesla & Schmittl(l2007l) . Most stars are Si- 
stars, but some of the hotter ones are also classified as He-weak, 
further all stars ha ve detected magnetic fields of a few hundred G 
up to several kG (Bvchkov et al. 2003). The use d stellar data is 
again taken from Kochukhov & Baanulo (2006), additional val- 
ues are from lLandstreet et al.l (120071) . In Fig. [8] we show a com- 
parison between the X-ray brightness and the stellar bolometric 
luminosity, errors are moderate, e.g. on log Lboi ~ 0.1. The com- 
parison sample reveals a striking similarity to our findings for 
IQ Aur and a 1 CVn; either the stars are X-ray bright with pre- 



dominantly log Lx > 29.0 erg s _1 or X-ray faint with most upper 
limits being around log Lx < 27.5 erg s , depending on the sen- 
sitivity of the respective exposure. Only high luminosity objects 
are detected in X-rays, supporting an intrinsic X-ray generating 
mechanism. Furthermore, all stars with a high luminosity are de- 
tected as X-ray sources without exception. Statistically the KS- 
probability of identical underlying distributions is well below 
1 % for our sample. A very similar dependence is present when 
using stellar mass, that is obviously related to luminosity; but 
notably, we find that strength of the magnetic field or the stellar 
age do not play a major role in making an Ap/Bp star an X-ray 
source (see Fig . [9j» . Unfortunately the sample does not cover the 
full parameter space, for example it lacks stars in the mass range 
3.3-3.8 M Q or with magnetic fields of B e ff > 5 kG and its size is 
still quite small. It thus does not allow to put strong constrains on 
the presence of a possible time-variability of the X-ray emission, 
but the fact that we detected all stars above a stellar luminosity of 
around 200 L (or 3.6 M ) and none below, suggests that a pos- 
sible transient aspect of the MCWS phenomenon might be of 
lesser importance. However, as is also obvious from Fig. [8] and 
Fig. [9] the X-ray luminosities of the detected stars have a large 
spread, indicating that other parameters play an additional role in 
determining the specific X-ray properties of the individual stars. 

5.4. Where does the X-ray emission of IQ Aur come from? 

The XMM-Newton observation clearly shows that the plasma 
temperatures of IQ Aur are significantly higher than those es- 
timated from the ROSAT data, where only weak indications for 
a hot component were found. This is very much relevant for the 
MCWS model, since Tx x V 2 h and thus high wind speeds and 
a strong magnetic confinement must be present. Temperatures 
of about 10 MK require typical shock velocities of several hun- 
dred up to 1000 km s . These values might be hard to obtain for 
AOp stars, though relatively high shock temperatures would be 
a natural explanation for the absence or weakness of a plasma 
excess at 2 MK. The relatively high O vm/O vn ratio in IQ Aur 
could be explained by a quite high temperature in the colliding 
wind shocks of Ap stars; for example at shock speeds of 500- 
600 km s" 1 noO vii excess emission is produced at all. At a few 
stellar radii a line driven wind has already nearly its terminal 
speed and for strong magnetic confinement a large longitudi- 
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nal velocity component is expected to be present, easily explain- 
ing the absence of a cool excess even in a wind shock scenario. 
Similarly, the X-ray luminosity of IQ Aur requires a rather high 
mass loss rate and a very efficient conversion of kinetic wind 
energy to X-ray emission. While not implausible, this requires 
IQ Aur to drive a strong, fast wind. However, given the presence 
of a disk, a fraction of the harder X-ray emission could be gener- 
ated in minor reconnection events. The O vn // /-ratio diagnostic 
indicates that the X-ray emitting plasma is not exposed to strong 
UV radiation. Thus, if not shielded by an unknown medium, the 
X-ray emitting regions are not predominantly close to the stellar 
surface, but rather at distances of several stellar radii. 

Our lower limit derived of from the X-ray spectra of d ~ 5 R* 
matches the value derived above for the inner disk radius, there- 
fore a disk-like confinement region at distances of 5 - 15 R* is 
a suitable candidate for the X-ray emitting region. On the other 
hand, the absence of significant rotational modulation over half 
of a rotation period implies, that the changing viewing geome- 
try introduces a negligible amount of time-variable, intervening 
material in the line of sight, caused either by an optically thick 
disk or the star itself. While the star as occulter appears unlikely 
given the large separation to the X-ray location, also for a disk 
rotational modulation is only expected if it is rather thick and 
the viewing angle changes significantly with phase. A small in- 
clination is unlikely, given the relatively large Vsini for an AOp 
star. However, if the magnetic obliquity is small, rotational mod- 
ulation effects are expected to be only minor. In this case the in- 
clination should be quite large and the system would be viewed 
nearly edge-on at all times. 

A strong flare-like event with 100 MK plasma is not ex- 
pected to be produced by wind shocks even for magnetically 
channelled winds, thus these X-rays must be generated by mag- 
netic reconnection. X-ray emitting structures on the surface of 
IQ Aur appear rather unlikely given the result from the X-ray 
diagnostics; additionally the magnetic field of Ap/Bp stars is 
overall a large scale one and structurally differs from coronae 
where most X-ray plasma is confined in small-scale structures. 
Therefore in the MCWS scenario the flare requires the presence 
of a rigidly rotating equat orial disk as pred icted by the above out- 
lined MHD simulations (ud-Doulaet a 1. 2008). In this scenario 
magnetic reconnection occurs, leading to the production of hot 
X-ray emitting plasma and further, episodic centrifugal break- 
outs of confined plasma produce flare-like events, originally pro- 
posed to explain X-ray flares on the Bp star <x Ori (see discus- 
sion below) and dud-Doula et alJ 120061) . In contrast to coronal 
flares, these magnetic reconnection events are driven by cen- 
trifugal mass ejections that disrupt outer disk regions. During 
breakout, the centrifugal forces on the accumulated disk mat- 
ter stretch the field lines more and more outwards and finally 
reconnection occurs. The event is accompanied by the ejection 
of outer disk mass, an inward propagating rebound and releases 
sufficient energy to produce hard X-ray flares. The finding of 
metallicity changes during the flare implies that the accumulated 
material in the circumstellar disk, which is heated during the dis- 
ruption, is more metal rich compared to the plasma that gener- 
ates the quasi-quiescent X-ray emission. Our rough estimations 
of the size of the flaring structure would imply that the outbreak- 
flare from IQ Aur was a more local event in the disk, involving 
a large clump or segment rather than the entire disk. Subsequent 
energy release or a change in the inner disk structure might be 
responsible for the flattening of the flare decay and the continu- 
ously enhanced level of X-ray emission after the event. Overall 
this scenario appears suitable to also explain the flare event on 



IQ Aur. In this case, we have likely observed the first X-ray flare 
from an AOp star. 



5.5. X-ray bright companions for Ap/Bp stars? 



As an alternative explanation it has been postulated that X- 
ray detected intermediate mass stars possess unknown late- 
type companions that actually produce the X-ray emission. 
Considering the estimated age of IQ Aur, i.e. several ten up to 
hundred Myr and comparing its X-ray luminos ity with studies 
of the 70 Myr old Pleiades (Mic ela et al.lll996l) . an active G or 
early-K dwarf as well as an M dwarf binary would be suitable 
candidate X-ray sources. However, young active stars have typi- 
cally significantl y higher plas ma temperatures than observed for 
IQ Aur (see e.g. Giidel 2004) and their X-ray light curves show 
quite frequent X-ray flares of factor up to a few, whereas IQ Aur 
is dominated by softer emission (Fig.[7J and exhibits except for 
the large flare only minor variability. Other properties like the 
high ///-ratio in the O vn triplet could be easily explained with 
a companion, typical densities are usually low and already a dis- 
tance of 0. 1 AU to IQ Aur would be in sufficient to completely 
erase the influence of the UV-field. Also the occurrence of a large 
flare are no problem for an active late-type companion and the 
associated abundance changes are reminiscent of the chromo- 
spheric evaporation scenario describing the flare evolution. 

Another criterion to consider is the fact that IQ Aur is also 
a bright radio source at 6cm with logL6 ~ 16 ergs -1 Hz -1 , 
whereas a 2 CVn with l og £6 < 14.9 ergs -1 Hz -1 is radio dark 
(radio properties from Linsk v et al.lll992l) . Since a rather tight 
correlation between X-ray and radio brightness measured at v > 
5 GHz is present for coronal sources, these measurements allow 
to test the companion scenario by checking the well established 
X-ray/Microwave(radio) correlation (also known as Giidel-Benz 
relation), which predicts Lx/Lr = l Q 15 - 5±0 - 5 Hz for magneti- 
cally active stars (iGiidel & Benj[l993l) . With a quasi-quiescent 
X-ray luminosity of logLx ~ 29.6 ergs -1 , IQ Aur is radio- 
overluminous by roughly two orders of magnitudes and thus 
significantly violates this correlation. The deviation is well in ex- 
cess of the typical observed scatter that is typically of a few. Note 
that similar X-ray and radio fluxes were repeatedly observed for 
IQ Aur, making intrinsic variability quite unlikely to be the rea- 
son for the discrepancy. Also the B8p He-w/Si star HD 144334 
(V929 Sco) with logL 6 * 16.0 ergs 1 Hz 1 and logL x ~ 
30.3 erg/s violates the Giidel-Benz relation by one to two or- 
ders of magnitude, similar to the B8V He-w star HD 146001 
with a radio detection at logLf, * 16.0 ergs -1 Hz -1 (plus one 
ILL. at logL6 ^ 15.9 ergs -1 Hz -1 ). An even stronger violation 
is observed for the B8IIIp He-w/Si star HD 142301 (V 927 Sco), 
with logL6 ~ 17.1 ergs -1 Hz -1 and logLx ~ 28.1 ergs -1 , this 
star is radio overluminous by more than four orders of magni- 
tude. On th e other hand, a ll non X-ray detected comparison stars 
studied bv iLinskv et all ([I992I i.e. HD 122532, HD 147010, 
remained also undetected at radio wavelengths at a limit of 
logLg < 15.7 ergs -1 Hz -1 . This points to a joint occurrence 
of strong X-ray and radio emission for the sample stars. The fact 
that IQ Aur and all other X-ray detected Ap/Bp stars violate the 
Giidel-Benz relation by orders of magnitudes suggests, that their 
emission does not originate from a late-type dwarfs, but a con- 
tributing companion in individual objects is not ruled out. 
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Fig. 10. Comparison of the quasi-quiescent X-ray spectra of the 
AOp star IQ Aur (black) and the B2p star cr OriE (red) using 
XMM/PN data. 



5.6. What about hotter magnetic stars? 

The MCWS model has also been applied to a variety of hot- 
ter stars and it is illustrative to compare them with IQ Aur. A 
well studied examp le is the magnetic 05 star l OriC, where 
iGagne et all d2005l) find that the dynamic MCWS model de- 
scribes the X-ray properties of the star very well. It has an 
average X-ray luminosity of Lx = 1 x 10 33 ergs -1 and uti- 
lizing among other diagnostics ///-line ratios, they conclude 
that the X-ray emitting plasma is located at distances of only 
0.2 . . .0.5 R* rather close to the surface of the star. l OriC ex- 
hibits rotational modulated X-ray emission, since it is a tilted 
magnetic rotator with / = 45° and /? = 42°, leading to a strong 
viewing angle dependent visibility and obscuration of the X-ray 
emitting torus. 9 l Ori C shows rather hard X-ray emission from 
very hot plasma at temperatures of up to 30 MK, but no strong 
flaring. Noteworthy, there have also been reports of X-ray flares 
on early type magn etic stars, a prominent exa mple is the B2p 
star cr OriE. While iGroote & Schmit tl (120041) attribute a flare 



seen with ROSAT to cr Ori E itself, iSanz-Forcada et al l (120041) 
report a long duration (> 10 h) flare with a count rate increase by 
a factor of ten in an XMM-Newton observation and concluded 
that it more likely originates from a companion. They derived a 
quasi-quiescent X-ray luminosity of Lx — 9 x 10 3l) ergs -1 and 
plasma temperatures up to 10-15 MK for cr OriE. Also non 
chemically peculiar, magnetic early-B stars are X- ray sources, 
e.g. th e B0.2 star t Sco with L x -2x 10 31 erg s' 1 dCohen et alj 
2003). The Chandra spectra show emission lines that are broader 
than for coronal sources, but surprisingly narrow for a stellar 
wind shock scenario; together with the presence of hard emis- 
sion and X-ray plasma locations at several stellar radii make 
also r Sco a candidate for MCWS generated X-ray emission. 
Anoth er example is the p ulsating, hot magnetic giant star Cep, 
where iFavata et al.1 (2009) again find a good agreement between 
the predictions from the dynamic MCWS model and X-ray ob- 
servations. It does again exhibit no flares and with temperatures 
well below 10MK the X-ray emission is rather soft, despite its 
X-ray brightness of log Lx — 4 x 10 30 erg s -1 . They derive dis- 
tances of m 5R* as the most likely location of the cooler X-ray 
emitting material, while the hotter plasma is likely closer to the 
star. Since /3 Cep has / = 45 and j3 = 85° but shows only minor 
modulation in its X-ray light curve, they exclude the presence 
of a thick disk and attribute the X-ray emission exclusively to 
magnetically channelled wind shocks. 

When comparing IQ Aur to the above mentioned massive 
magnetic stars, we find some similarities and overall trends. 



Clearly, in general the X-ray luminosity decreases towards less 
massive stars, however all magnetic stars are with log Lx/^boi ~ 
-6 X-ray overluminous compared to similar 'normal' massive 
stars. The bulk of the X-ray emitting plasma is located further 
away from the central star going from O over B to A-type stars 
and likely in the direction of increasing 77*. IQ Aur is also by far 
the coolest object and the only star in this comparison, where the 
// / ratio of O vn is high. Comparing the X-ray temperatures, we 
find that the plasma on IQ Aur is cooler than on 9 l Ori C, compa- 
rable to cr Ori E, but hotter than on fi Cep. For stars with a mod- 
erate wind speeds well below 1000 kms -1 , in the MCWS sce- 
nario the lower temperature plasma (< 5 MK) is predominantly 
produced by wind shocks, whereas hotter plasma (> 10 MK) is 
generated by reconnection in the disk. Notably, the strong flar- 
ing sources IQ Aur and cr Ori E have by far the largest confine- 
ment parameter (log r] m ~ 6 ... 7), favoring the effective accumu- 
lation of matter in the equatorial region and consequently break- 
out events. As shown in Fig. [10] also their quasi-quiescent spec- 
tra are very similar, despite the fact that underlying stars differ 
significantly and cr Ori E has a twenty times higher X-ray lumi- 
nosity. 

5. 7. The Ap/Bp star X-ray puzzle 

Our findings provide some clear hints for the solution to the 
X-ray puzzle of AOp stars. We do not find strong support for 
the hypothesis that the X-ray emission of Ap/Bp stars is pro- 
duced by X-ray emitting companions, but cannot rule out that 
in some cases multiple sources contribute to the observed emis- 
sion. It is much more likely that the X-rays originate from intrin- 
sic emission mainly generated in magnetically channelled wind 
shocks. However predominantly or even exclusively the more 
massive, high luminosity stars appear as X-ray sources. These 
late B/early A stars drive a stellar wind with a sufficiently large 
mass loss rate, whereas further stellar properties like the mag- 
netic topology or orientation and transient phenomena influence 
primarily the details of their X-ray appearance. In the case of 
very strong magnetic confinement an equatorial disk structure 
forms, where magnetic reconnection and centrifugal breakout 
events occur that produce hard X-ray emission and occasional 
flare-like events. 

Dedicated modeling of AOp stars and their magnetically 
channelled winds, companion searches, as well as more exten- 
sive studies of similar objects are highly desired to gain even 
deeper insights in the origin of X-ray emission in these remark- 
able stars. 

6. Summary and conclusions 

1 . We detect X-ray emission from IQ Aur at log Lx = 
29.4 ergs -1 , but failed to detect a 2 CVn at logL x < 
26 erg s -1 although both are AOp stars. The quasi-quiescent 
spectrum of IQ Aur exhibits significant emission from in- 
termediate temperature and even hot plasma (> 10 MK), in 
addition to a strong cool component at 3 MK. The UV sen- 
sitive Ovu ///-ratio is high, thus the X-ray emission very 
likely originates from regions at distances of > 5R* above the 
stellar surface. The X-ray light curve shows minor variability 
without clear signs of a viewing angle dependent rotational 
modulation. 

2. A massive X-ray flare with rise time of about 1 ks and fast 
decay occurred during the IQ Aur observation. The X-ray lu- 
minosity at the flare peak is about logLx = 31.5 ergs -1 and 
plasma temperatures approach 100 MK. The flare originates 
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in a rather compact region and increases the metallicity of 
the X-ray emitting plasma significantly. 

3. Our X-ray data of IQ Aur can be reconciled within a dy- 
namic MCWS model with strong magnetic confinement and 
build up of a rigid disk where magnetic reconnection events 
occur. In this scenario, the flare-like events was produced by 
the centrifugal breakout of confined material. Alternatively, 
IQ Aur might be a multiple system where likely both com- 
ponents contribute to the observed X-ray emission. 

4. Magnetically channelled stellar winds are the most promis- 
ing explanation for the X-ray detection of Ap/Bp stars, al- 
though the presence of X-ray emission strongly depends on 
stellar parameters. A suitable quantity for setting the stage of 
MCWS generated X-rays appears to be a high stellar bolo- 
metric luminosity that supports a high mass loss, but mag- 
netic topology, orientation or transient phenomena also play 
a role in determining the specific X-ray properties of an in- 
dividual Ap/Bp star. 
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